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The nonlinear response to an external electric field is studied for classical non-interacting charged 
particles under the influence of a uniform magnetic field, a periodic potential, and an effective 
friction force. We find numerical and analytical evidence that the ratio of transversal to longitudinal 
resistance forms a Devil's staircase. The staircase is attributed to the dynamical phenomenon of 
mode-locking. 



The electron transport in two-dimensional periodic 
arrays of scatterers has been actively studied for the 
last decade. One of the most interesting feature is the 
plateau-like behaviour in Hall resistance as well as peaks 
in magnetoresistance [Q,^ at low magnetic fields below 
the quantum Hall regime. The peak structure in magne- 
toresistance has been attributed to the electron cyclotron 
orbits which enclose an integer number of scatterers . 
The cyclotron motion is important when the electron 
mean free path I {I is measured in the absence of peri- 
odic scatterers) is greater than the period a of the regular 
scatterers. 

In this Letter, we present a theory for electron trans- 
port in the deeply diffusive regime I <gC a where cyclotron 
motion is not relevant. We predict a new and interesting 
effect where the structure in the magnetoresistance is as- 
sociated with fractional numbers. We will show that the 
relative ratio between the transversal and longitudinal 
magnetoresistance forms a Devil's staircase. These frac- 
tal staircases originate from a dynamical phenomenon 
known as mode-locking which appears naturally in the 
context of circle maps ^ . Our prediction is based on the 
finding that the dynamics in a two-dimensional periodic 
array of scatterers with intrinsic momentum relaxation 
effectively reduces to a circle map. 

Let us consider first a particle of mass m, charge 
q in crossed electric {E = Ex) and magnetic fields 
(B = —Bz) with intrinsic momentum relaxation. The 
resulting drift motion x{t) = (a;(t), y{t)) can be described 
by adding a frictional force, proportional to the drift ve- 
locity V, to the equations of motion yielding 



X 

y 



mix = —qBvy + qE — mVx/T 
mVy — qBvx — mVy/r 



(1) 



with T being the momentum relaxation time. Cyclotron 
motion only exists for a transient time; each trajec- 
tory finally converges to a straight line xit — > oo) = 
a;(0) -I- {v)t, with the time-averaged velocity {v) = 
(m/r, qB)qE/ {q^B^-\~m'^ /t'^). With the charge density n 
and the current density J — qn{v), it is easy to calculate 
the resistivities defined hy E ~ pj and Onsager's rela- 
tions pxx — Pyy and Pyx — —pxy- the diagonal resistivity 
Pxx = mjq^nT is constant, whereas the off-diagonal re- 
sistivity pxy = B/qn depends linearly on the magnetic 



field - no plateaus are present. 

We now add a periodic potential V = VqF{x) with \F\ 
of order one to the equations of motion (^, assuming 
that the periodicity of the potential, a, is much larger 
than the length scale related to the intrinsic momentum 
relaxation. When we scale the coordinates x with a and 
the time t with r, we get the same equations of motion 
for a given set of three dimensionless parameters: B = 
qTB/m, E = qT'^E/{ma), and Vq = Vo/{qEa). We use 
the rescaled magnetic field B as a variable parameter and 
fix the rescaled electric field and potential strength to 
E = 15.3 and Vq = 0.12 if not otherwise stated. Though 
we could choose various sets of parameters (r, a, E, etc.) 
for the chosen values of E and Vq , for definiteness we will 
take T = 7.6 10-"^^ e = -2.3 W^V/cm, and a = O.S^im, 
and fix these for the rest of the paper. The particle is an 
electron in GaAs sample with q = — e and to = 0.067TOe. 
By taking a typical value of Fermi velocity « 3 x 
10^cto/s, the mean free path of our system is I — VfT w 
0.02^m < a. 

First, we show how the particle's motion is related to a 
one-dimensional circle map through analytical consider- 
ations for sufficiently small potential strength. Note that 
our numerics do not rely on this requirement. Due to the 
periodicity of the potential it is sufficient to consider the 
unit cell of the potential with periodic boundary condi- 
tions, i.e. we can compactify the (x, ?/)-plane to a two- 
dimensional torus. A typical asymptotic solution for zero 
potential strength as discussed above appears here as a 
quasiperiodic trajectory filling densely a two-dimensional 
invariant torus v = constant in four-dimensional phase 
space. Trajectories starting away from the torus are 
attracted towards it and asymptotically converge onto 
it. For finite potential strength the torus is smoothly 
deformed or collapsed into lower-dimensional objects or 
broadened to a higher-dimensional object (loosely speak- 
ing, a torus with finite but very small thickness). The 
latter possibility can be ignored to a very good approx- 
imation as we will see later. A point on the torus (and 
on the lower-dimensional objects) is uniquely labelled 
by unit-cell coordinates (xi mod 1, 2:2 mod 1) defined by 
X = xiai -\- X2a2, with ai = (an, 012), 02 = (021,022) 
being the lattice vectors of the unit cell. The long- 
time behavior is therefore completely described by two- 
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dimensional dynamics [xi{t),X2{t)]. As usual, such dy- 
namics can be conveniently investigated by introducing 
Poincare surfaces of section defined by xi{t„) mod 1 = 0, 
where n is an integer and modulo 1 restricts the variable 
to the interval [0, 1). 

The discrete evolution of X2.n = X2{tn) is then gov- 
erned by a one-dimensional map X2,n+i — f{x2,n) where 
/ is a smooth and monotonic function. Due to the pe- 
riodicity of the potential the condition / (0:2) modi = 
f{x2 mod 1) holds. Such a map is called an invertihle 
circle map. Its rotation number, 



p = lim 



X2,n - X2.0 



n 



(2) 



is well-defined and independent of the initial condition 
2^2,0 @1- A rational value oi p = p/q {p and q are inte- 
gers without common divisor) indicates periodic motion, 
whereas an irrational value indicates quasiperiodic mo- 
tion; chaotic motion is not possible. 

We now derive a simplified circle map valid in the limit 
Vq ^ I. If Vo = then the asymptotic dynamics obeys 
Vx = i'y — 0, leading to f{x2) — X2+P0 with the constant 
Pq{B); here we have simply p — po, so po can be regarded 
as the unperturbed rotation number. The asymptotic 
dynamics for nonzero but sufficiently small Vq is over- 
damped, i.e. \vx\ ^ b^l/''' and l^j^l <C \vy\/T. As in the 
potential-free case the number of differential equations 
reduces to two, enabling us to determine / up to first 
order in Vq 



f{x2,7i;Po,E,VQ) « X2,n + Pa + Vo g{x2,n\Pa) (3) 



with 



/Oo("22 + ^21) + 2/Oo(aiia2i + 012022) 



all 
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X2=P0Xl+X2,: 



(4) 



5 is a nontrivial periodic function of X2,n- Note that / is 
independent of the parameter E in the present approx- 
imation. Using Stoke's theorem it can be shown that 

gdx2 — 0. For circle maps of the form given in Eq. (j^) 
with the above-mentioned properties of g it is proven that 
for each Vq there exists a Devil's staircase, a monotoni- 
cally increasing function p — p{po) with plateaus around 
each rational value of po. The width of the plateaus, 
Apo, is proportional to Vb @- Similarly, the widths of 
the plateaus, AB, of the function p — p{B), are propor- 
tional to Vb since AB « {dB/dpo)Apo. 

The rotation number p is directly related to the 
continuous-time dynamics. This can be seen by rewriting 
Eq. (0) with V — viai + V2a2 as 



(^^2) _ -012 + aii{vy)/{vx) _ -012 -I- aii(Jxy/crx 

(Vi) 022 - a2l{Vy)/{Vx) a22 ~ a2iaxy/(7xa 



Hence, p is given by the ratio of the conductivities. In the 
case that the periodic potential is invariant under n/m- 
rotations, to = 2, 3, . . ., Onsager's relations show that p 
can also be expressed by the ratio Pxy/Pxx- 

We now present numerical evidence showing that for 
finite Vb < 0.4 the continuous-time dynamics is indeed 
given by a one-dimensional circle map and that the mag- 
netoresistance shows a fractal structure. As an example 
of the periodic scatterers, we take a potential with hexag- 
onal symmetry shown in Fig. |^a, 

V{x, y) = Vb I cos (^(\/3x - y)) cos y 

cos (^^{VSx + y)^ cos (■^(\/3: 

(liVSx + y) 
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cos I — y ) cos 



x-y) 
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with /3 = 8. The lattice vectors are ai — {\/3, 0) and 
12 = (— \/3/2, — 3/2). We have used the Runge-Kutta 
method |^ with 10^ time steps of variable size to inte- 
grate the equations of motion. We found that all initial 
conditions, for fixed parameters, lead to the same mean 
velocity {v). So only one single orbit is needed to cal- 
culate the current density and the resistivities without 
employing the Kubo formula [Q. This already indicates 
the collapse to an invertible circle map. To show the re- 
duction to the map explicitly, we calculate a sequence of 
points {x2,n mod 1, a;2,n+i mod 1) by solving the complete 
set of differential equations for different initial velocities 
from one Poincare section to the next. Figures ^ and 
b show that these points indeed lie very close to a line, 
so a;2.n+i is to excellent approximation a function of X2^n 
independent of the initial velocity. We have computed 
{x2,n+i—X2,n — Po)/Vo for different parameters Vb and E. 
It can be seen from Fig. [^c that this quantity is roughly 
independent of Vb, even though small deviations indicate 
weak nonlinear behavior in Vb- Figure |2|d shows that 
{x2,n+i — X2.n — Po) /Vo docs not chaugc when E is var- 
ied over two orders of magnitude. Figures |^-d therefore 
confirm the reduction to the circle map and its parameter 
dependence as predicted in Eq. (^. 
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FIG. 1. a) Contourplot of the potential in Eq. (|6|) with /3 : 
b) Orbit with B = 8.7T and p = 1 in the unit cell, c) B = 9.77T 
and p = 1. d) 9.8T and p = 1.01407.... e) 2.9T and p = 1/2. f) 
4.4r and p = 2/3. 
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FIG. 2. a) X2,„+imodl vs. X2,n for B = 8.7T. b) 9.77T. c) 
{X2,n+l — X2,n — Po)/Vb VS. X2,n for _B = 8.7T, Vb = 0.12 and 
E = 15.3 (black); Vq = 0.06 (red); Vb = 0.03 (green), d) % = 0.12 
and E = 15.3 (black); E = 1.53 (red); E = 0.153 (green). 

Figure H is our main result. It shows the Devil's stair- 
case p vs. pq. Both quantities are computed from the 
resistivities via the relations p = 2/(1 + V^pxx/ Pxy) and 
Pa = 2/(1 - VS/B) from Eq. (|) with the lattice vec- 
tors of the rotational symmetric potential (^). Note the 
exactness of the plateaus. Let us discuss some features 
of the Devil's staircase in relation to the continuous and 
discrete-time dynamics. Figures |^b and c show two orbits 
moving in slightly different magnetic fields. Both have 
synchronized velocities (vi) = (^2), i.e. p = 1. Each orbit 
with p = 1 is a stable fixed point = f{xl) in the corre- 
sponding map (red circles in Figs. ^ and b; blue circles 
mark unstable fixed points). This synchronization phe- 
nomenon, usually called mode-locking in nonlinear dy- 
namics, is therefore due to the robustness of fixed points 
under variation of a parameter, here po{B). Orbits with 
p ^ 1 are in a different mode; in particular, orbits with 
irrational rotation number are not periodic. In Fig. 
we see an orbit with irrational p close to one. It stays 
for a long time near a periodic orbit with p — I but from 
time to time it escapes, thereby filling the entire unit cell 
densely. This phenomenon is called intermittency In 
terms of the circle map, the transition from p = 1 to irra- 
tional p is a saddlenode bifurcation: as B is varied, stable 
and unstable fixed points come closer and closer as illus- 
trated in Figs. ||a and b, and finally destroy each other, 
yielding quasiperiodic motion. Varying B further can 
lead again to periodic motion. Figures ^ and f give two 
examples with p = 1/2 and p = 2/3. The corresponding 
map has a stable fixed point of higher period, e.g. for 
p = 1/2 the fixed point condition is X2 — /[/(a;^)]. 
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FIG. 3. Rotation number p vs. po- Arrows assign the rational 
value of the rotation number p = p/q to each plateau. The inset 
highlights the fractal structure of the Devil's staircase. 

We also observe a Devil's staircase- like function pxy{B) 
in Fig. W, which reminds one of the quantum Hall ef- 
fects [plJlOl], but here the 'plateaus' are not perfectly flat. 
One may note also that the way p^x varies for a num- 
ber of p^, J,- 'plateaus' has similarities to the quantum Hall 
effects: going from the centre of a pa;j,-'plateau' towards 
its border, pxx increases and reaches its maximum ex- 
actly when the 'plateau' ends. This can be explained 
qualitatively. First, note that an orbit from the centre of 
the p = 1-plateau (Fig. |l|b) avoids the potential's steep 
maxima at the corners of the unit cell. Hence, the di- 
agonal resistivity roughly equals the potential-free value 
Pxx = m/q^TiT « 3.12fcr2. We infer from Figs. |^b and 
c that the closer B is to the border of the plateau the 
closer is the orbit to the potential's maxima. This effec- 
tively slows down the orbit, which explains the increase 
of Pxx (and Pxy)- Leaving the plateau finally reduces the 
resistivities because the intermittency effect becomes less 
pronounced with increasing |p — 1|; see Fig. ^d. 
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FIG. 4. OfT-diagonal resistivity pxy and diagonal resistivity pxx 
vs. magnetic field B. Arrows assign the rational value of the rota- 
tion number p = p/q to each plateau. 

The example in Fig. |^ confirms that the widths of the 
plateaus, AB, is roughly proportional to Vq as predicted 
from the analysis of the simplified circle map in Eq. (|^). 
By noting that denotes the dimensionless potential 
strength we see an interesting analogy with the quantum 
Hall effect where the size of plateaus is sensitive to the 
disorder strength of the systems [|n| . Here we have the 
advantage of the use of the simplified circle map to anal- 
yse the plateau size, although the overall structure of the 
hierarchy of the plateaus is also too complicated to be 
dealt with in an analytic fashion. 



pies for fixed applied voltages, one needs to use metallic 
leads which span the entire length of two opposite edges 
of the sample [ p^ . 

In summary, we have calculated the magnetoresistance 
of a lateral surface superlattice with strong momentum 
relaxation where the cyclotron motion is not involved. 
Our calculations show a fractal plateau structure in the 
magnetoresistance which stems from purely classical non- 
linear dynamics. We have explained our calculational 
results in terms of the theory of circle maps where the 
Devil's staircase is already well understood. 

We acknowledge helpful discussions with P. Fulde, 
K. V. Klitzing, R. Gerhardts, K. Richter and R. Klages. 
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FIG. 5. Dependence of the p = 1-plateau on Vb; cf. Fig. ^. 

Since the particle can be trapped in the very flat local 
minima of the potential (H), a finite dc electric field is 
necessary to overcome the potential trap as in the ex- 
periment in Ref. ||l^. The threshold |13| value Et of 
the electric field for the finite electrical current might be 
crudely estimated by Et ^ Vq/ (qa) , which is in agree- 
ment with our calculations (not shown). However, in 
reality the threshold electric field will be much smaller in 
the presence of the degenerate electron gas because the 
weak local trap potential will be screened [Q. While 
we have chosen in our numerical calculation a large elec- 
tric field (2.3 10''V/cm) and a low mobility of the sample 
{fi = er/m « 20000cm^ /Vs), the Devil's staircase may 
be observed in a broad range of parameters. We also note 
that the specific form of the potential (^) is not relevant. 
However, the steepness of the potential influences the 
shape of the plateaus. Large steepness (large /3) results 
in wide and curved plateaus. It should be mentioned also 
that the Devil's staircase will not be seen in the conven- 
tional Hall bar where the electrical current is fixed along 
the bar and the induced Hall voltage is measured. To 
measure the electrical current in two-dimensional sam- 
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